Black hole state counting in Loop Quantum Gravity: A number theoretical approach 
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We give an efficient method, combining number theoretic and combinatorial ideas, to exactly 
compute black hole entropy in the framework of Loop Quantum Gravity. Along the way we provide 
a complete characterization of the relevant sector of the spectrum of the area operator, including 
degeneracies, and explicitly determine the number of solutions to the projection constraint. We use 
a computer implementation of the proposed algorithm to confirm and extend previous results on 
the detailed structure of the black hole degeneracy spectrum. 
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Any proposed quantum theory of gravity must account 
for the states responsible for black hole entropy. Within 
Loop Quantum Gravity (LQG) entropy can be studied by 
using the Isolated Horizon framework, The counting 
of states is reduced in this setting to a well defined com- 
binatorial problem. It gives rise, in the asymptotic limit, 
to the semiclassical Bekenstein-Hawking formula 0, Q 
corrected to the next relevant order by a term logarith- 
mic in the area. A computer assisted study has been 
carried out for small black holes up to two hundreds of 
Planck areas £p Q. This has unearthed a very inter- 
esting behavior in their degeneracy spectrum, namely, 
an equidistant "band structure" with important physical 
consequences. The most relevant of them is the effec- 
tive quantization of black hole entropy ^ . A qualitative 
understanding of the origin of this behavior has been ob- 
tained in . However, no detailed theoretical description 
of this phenomenon has been available to date owing to 
the incomplete characterization of the area spectrum, on 
one hand, and the lack of exact manageable solutions 
for some combinatorial problems (involving the so called 
projection constraint), on the other. 



standing of the different factors that come into play to 
reproduce the features previously observed in the black 
hole degeneracy spectrum. In addition, it can be effi- 
ciently used to perform exact entropy computations - 
extensible to large areas- that improve and confirm the 
results obtained by brute force methods in Q. 

We start by characterizing the area eigenvalues and 
their degeneracies. In LQG the black hole area is given 
by an eigenvalue A of the area operator 
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In this letter we present a satisfactory solution to both 
types of difficulties, giving a precise characterization of 
the area spectrum by relying on number-theoretic meth- 
ods, and addressing the combinatorial problems related 
to the projection constraint. We do it for the original 
counting of states proposed in [3 and carried out in , 
and also for the one described in p] . The method that we 
discuss in the following will allow us to have a full undcr- 



where 7 denotes the Immirzi parameter. Notice that 
these do not give the full area spectrum but, for the case 
of isolated horizons relevant here we only need equation 
([1]). The labels ji are half-integers, ji G N/2, associ- 
ated to the edges of a given spin network state. They 
pierce the horizon at a finite set of N distinguishable 
points called punctures Horizon quantum states are 
further characterized by an additional label mj. In the 
case where we have spherical symmetry a projection con- 
straint 
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must be satisfied by the mj. There are two inequivalent 
proposals in the literature to account for the relevant 
microscopic configurations 0, Q • When taken as a purely 
combinatorial problem, they differ in the range of the 
label m/. In the standard (DLM) counting performed 
in Q one takes to/ G {~ji,ji}, whereas the counting 
proposed in Q (that we will refer to as the GM counting) 
assumes that to/ can take all the allowed values for a spin 
component to/ g {-j/, -j/ + 1, . . . , j/ - 1, j/}. 

A first problem that we address is the characterization 
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of the numbers belonging to the spectrum of the area op- 
erator restricted to the vector subspace spanned by spin 
network states having no vertices nor edges lying on the 
black hole horizon. In the following when we talk about 
the area spectrum wc refer, in fact, to this restriction. 
The first question that we want to consider is: Given 
A S M, when does it belong to the spectrum of the area? 
In order to simplify the algebra and work with integer 
numbers we will write ji = ki /2 in the following, so that 
the area eigenvalues become 

1=1 k=l 

Here we have chosen units such that Anjij^ — 1, and 
the rifc (satisfying ni + ■ ■ ■ Uk-^^^ = N) denote the num- 
ber of punctures corresponding to edges carrying spin 
k/2. An elementary but useful comment is that we can 
always write yj [k + 1)^ — 1 as the product of an inte- 
ger and the square root of a square-free positive integer 
number (SRSFN) by using its prime factor decomposi- 
tion. Hence, with our choice of units, only integer linear 
combinations of SRSFN's can appear in the area spec- 
trum. The questions now are: First, given such a linear 
combination, when does it correspond to an eigenvalue 
of the area operator? If the answer is in the affirmative, 
what are the permissible choices of k and compatible 
with this value for the area? 

In the following wc will take advantage of the fact that 
SRSFN's are linearly independent over the rational num- 
bers (and, hence, over the integers) i.e. qi^/pi + • • • + 
<lr\fVr — 0, with G Q and different square- free in- 
tegers, implies that qi — ^ for every i = 1, . . . ,r. This 
can be easily checked for concrete choices of the and 
can be proved in general (see for instance [^). We can 
answer the two questions posed above in the following 
way. Given an integer linear combination of SRSFN's 
Si=i liy/Wi^ where € N, we need to determine the val- 
ues of the k and rik , if any, that solve the equation 

£ nfe\/(fc + 1)2 - 1 = ^ q.^^. (3) 

fe=l i=\ 

Each •\/(fc + 1)2 — 1 can be written as an integer times a 
SRSFN so the left hand side of ^ will also be a linear 
combination of SRSFN with coefficients given by integer 
linear combinations of the unknowns n^. As a prelimi- 
nary step, let us find out -for a given square-free positive 
integer p^- the values of k satisfying 

V(fc + l)2-l = yVP^, (4) 

for some positive integer y. This is equivalent to solving 
the Pell equation — piy^ = 1 where the unknowns are 
X := fc + 1 and y. Equation ((4]) admits an infinite number 
of solutions {kl^,yl^), where m € N (see, for instance, 
n). These can be obtained from the fundamental one 



{k\, yl) corresponding to the minimum, non-trivial, value 
of both kl^ and y^. They are given by the formula 

fcm + 1 + y'mVp^ =ik\ + i + 2/ivp^)™. 

The fundamental solution can be obtained by using con- 
tinued fractions ^2]. Tables of the fundamental solution 
for the smallest pi can be found in standard references 
on number theory. As we can see both /c^ and grow 
exponentially in to. By solving the Pell equation for all 
the different pi we can rewrite ([3]) as 

r oo r 
i—l m—1 i—l 

Using the linear independence of the the previous 
equation can be spht into r different equations of the type 

oo 

y'ni'^K. = * = (5) 

m—1 

Several comments are in order now. First, these are dio- 
phantine linear equations in the unknowns rij,i with the 
solutions restricted to take non-negative values. They 
can be solved by standard algorithms (for example the 
Frobenius method or techniques based on the use of 
Smith canonical forms). These are implemented in com- 
mercial symbolic computing packages. Second, although 
we have extended the sum in ([5]) to infinity it is actu- 
ally finite because the grow with m without bound. 
Third, for different values of i the equations (O are writ- 
ten in terms of disjoint sets of unknowns. This means 
that they can be solved independently of each other -a 
very convenient fact when performing actual computa- 
tions. Indeed, if {k^rnxiVmx) (^ma'J^ma) ^''^ solutious 
to the Pell equations associated to different square-free 
integers andpi2; then k^^^ and A;^^ must be different. 
This can be easily proved by reductio ad absurdum. 

It may happen that some of the equations in ([5]) ad- 
mit no solutions. In this case J2i=i 1i\/Pi does not be- 
long to the relevant part of the area spectrum. On the 
other hand, if these equations do admit solutions, the 
^i\/Pi belong to the spectrum of the area operator, 
the numbers k\^ tell us the spins involved, and the n^i 
count the number of times that the edges labeled by the 
spin A:*„/2 pierce the horizon. A set of pairs {(fc^,nj,i )} 
obtained from the solutions to equations ([3]), (|4]), and ([5]) 
will define what we call a spin configuration. The number 
of different quantum states associated to each of these is 
given by two degeneracy factors, namely, the one coming 
from reorderings of the fc/-labels over the distinguish- 
able punctures (r-degeneracy) and the other originating 
in all the different choices of TO/-labels satisfying (2), (m- 
degeneracy) . The combinatorial factors associated to the 
r-degeneracy are straightforward to obtain and appear in 
the relevant literature. 

Let us consider then the m-degeneracy. The problem 
that we have to solve is: Given a set of (possibly equal) 
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spin labels ji, I = 1, . . . , iV, what are the different choices 
for the allowed mj such that ^ is satisfied? Notice that 
an obvious necessary condition for the existence of solu- 
tions is that ji £ I^- 

In the standard DLM approach the number of differ- 
ent solutions for the projection constraint can be found 
by solving the following combinatorial problem (closely 
related to the so called partition problem): Given a set 
IC — {ki, . . . , fcjv} of iV -possibly equal- natural numbers, 
how many different partitions of /C into two disjoint sets 



/Ci and IC2 such that ^ 



keKi 
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k do exist? The 



answer to this question can be found in the literature 
(see, for example, Q and references therein) and is the 
following 
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(6) 
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where M = 1 + E/=i ''■J- This expression can be seen 
to be zero if there are no solutions to the projection con- 
straint. 

Let us consider now the GM proposal. The problem 
is equivalent in this case to counting the number of ir- 
reducible representations, taking into account multiplic- 
ities, that appear in the tensor product ^f^iiji], where 
[ji] = denotes the irreducible representation of 

SU{2) corresponding to spin jj. In order to solve this 
problem we rely on techniques developed in the context 
of conformal field theories (see al so [To| ) and in the 
spectral theory of Toeplitz matrices The starting 

point is to write the tensor product of two SU(2) repre- 
sentations in the form 
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where the integers J^kik2 > called fusion numbers , tell 
us the number of times that the representation labeled by 
/c3/2 appears in the tensor product of [A:i/2] and [^2/2]. 
For each k G NU{0}, we introduce now the infinity fusion 

matrices {Cu)kik2 '■— -^kik' where fci, ^2 € NU{0}. These 
can be shown to satisfy the following recursion relation 



fe+2 
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Ck 



0,1,. 



(7) 



where we have introduced the notation X := Ci. Explic- 
itly Xkj^k2 = Ski,k27zX,+ ^fci,fc2+i, which shows that X is 
a Toeplitz matrix [11]. The solution to ([7]), with initial 
conditions Cq ~ I and Ci — X, can be written as 

Ck = Uk{X/2), A: = 0,1,... 

in terms of the Chebyshev polynomials of the second kind 
Uk- The tensor product of an arbitrary number of rep- 
resentations can be decomposed as a direct sum of ir- 
reducible representations by multiplying the fusion ma- 
trices introduced above. By proceeding in this way we 
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^(Cfe^Cfeg • • •Cfc„)feifc 
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Notice that the product of matrices appearing in the pre- 
vious formula is, in fact, a polynomial in X. The total 
number of representations, that gives the solution to the 
combinatorial problem at hand, is simply given by 



^(Cfc2Cfe3 • • • Cfe„) 



kik • 



(8) 
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This is just the sum of the (finite number of non zero) el- 
ements in the fci row of the matrix Ck2Ck^ ■ ■ ■ Cfe„. A 
useful integral representation for this sum can be ob- 
tained by introducing a resolution of the identity for 
X as in [ll| and the well-known identity C/„(cos6') — 
sin[(n -I- 1)0]/ sin 6*. In fact, equation ([8]) can be equiva- 
lently written as 



— i AO cos - 







cos - — COS (yK+ -) 
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sin(fc/ + 1)9 



, (9) 



where K = ki + ■ ■ ■ + k^- This is related to the well 
known Verlindc formula for SU{2) 

The procedure to calculate the black hole spectrum 
described above can be efficiently implemented in a com- 
puter, for instance using Mathematica. This allows us 
to analyze in detail the different factors that shape the 
degeneracy spectrum. First of all, the fact that the dio- 
phantine equations are decoupled allows us to obtain 
the configurations compatible with a given value of area 
A = El=i Qi\/Pi the cartesian product of the sets of 
solutions to the diophantine equations for each pi . Let us 
then begin by analyzing the results for area values of the 
form A = q^, with g e N and a fixed SRSFN. What 
we see in this case is that the r-degeneracy -coming from 
the reordering of puncture labels- will be maximized by 
those configurations having both a large number of dif- 
ferent values of k and a large number of punctures. For 
a fixed area value these two factors compete with each 
other because higher values of k imply a lower number of 
punctures. On the other hand the m-degeneracy shows 
an exponential growth with area (both in the DLM and 
GM countings). When the two sources of degeneracy are 
taken into account -in the present case involving a single 
SRSFN- the total degeneracy can be seen to be dom- 
inated by the m-degeneracy. The reason for this domi- 
nance of the m-degeneracy is that the number of different 
(small) values of k available within the set of solutions to 
the Pell equation for a given p is limited, and hence only 
a few possibilities of reordering exist. 

This situation is expected to change drastically when 
we consider areas A = E[=i Qi\/Pij with r > 1, built as 
linear combinations of different SRSFN's. In this case 
it is possible to obtain configurations with a large num- 
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ber of different small values of k (associated to differ- 
ent SRSFN's). The effect of considering linear combi- 
nations involving several SRSFN's produces a very dis- 
tinctive feature when the r-degeneracy is plotted as a 
function of area, namely, it creates a "band structure" 
where high values of degeneracy alternate with much 
lower ones. Furthermore, maxima and minima are evenly 
spaced. When this behavior is considered together with 
the m-degeneracy we obtain the regular pattern shown 
in Figure [TJ 
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FIG. 1; Plot of the black hole degeneracy as a function of 
the area (expressed in units of to facilitate the comparison 
with the results obtained in Q]). 



Several remarks are in order now. First, we want to 
point out that the result obtained from the explicit com- 
putational analysis carried out in [1] (by using the GM 
counting) is exactly recovered with the new approach. 
The fact that the same result is obtained from two com- 
pletely independent procedures (a brute force approach 
and the algorithm proposed here) provides strong evi- 
dence for the reliability of both computations. Second, 
the structure of the degeneracy spectrum obtained by 
using the DLM and GM countings is basically equal. 
They differ only in the absolute values of the degener- 
acy whereas the band structure (including the position 
and spacing of the bands) is the same. This can be un- 
derstood in our framework because the terms account- 
ing for the r-degeneracy, responsible for this effect, co- 
incide for both counting procedures. This justifies the 
appearance of the constant x obtained in 0, @. Third, 



once we understand how the r-degeneracy works, we see 
that the area values for which the degeneracy is large are 
those that can be written as linear combinations of the 
SRSFN's originating from small solutions k to the cor- 
responding Pell equation. Thus, considering these linear 
combinations will suffice to account for the band struc- 
ture. The remaining area values give rise only to very 
low degeneracies. 

Summarizing, we have been able to find a number- 
theoretic/combinatorial way to tackle the problem of 
calculating the degeneracy spectrum of spherical black 
holes in LQG. Our procedure has several advantages 
over previous approaches. First, we have been able to 
characterize the area spectrum in a proper way, giving 
an algorithm to explicitly find every single spin con- 
figuration contributing to each value of the area spec- 
trum. In particular, the degeneracies of the area eigen- 
values can be obtained. This has allowed us to repro- 
duce and understand the band structure already ob- 
served in [3] for the black hole degeneracy spectrum in 
a much more efficient way. We not only recover previ- 
ous results obtained by using a brute force algorithm, 
but easily extend them to area values significantly larger 
than those reached in [3| (see Figure [T]). Moreover, 
with our methods it is possible to compute the config- 
urations and degeneracy even for much larger values of 
area. As a token we give the degeneracy for an area 
of 8320^/2 + 14400^3 + 2240^6 + 4640a/15 + 1120^35, 
which is 3.46437296507975 • • - x lO^^'^ao, Finally, the con- 
crete procedures and explicit formulas given in the letter 
offer a good starting point to study the asymptotic be- 
havior of the entropy as a function of the area of a black 
hole. This could help us investigate whether the effective 
entropy quantization discussed here is present in macro- 
scopic black holes. 
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